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Abstract. We model the positional dependence of the 
optical polarization (BVRI-bands) in the circumstellar 
disc of (3 Pictoris as observed by Gledhill et al. (1991) 
and Wolstencroft et al. ( 1995| ). The particles are spheri- 



cal, have a size distribution n(a) ~ a~ q and their num- 
ber density decreases with distance from the star as r~ s . 
We consider both compact and porous grains. Varying the 
grain size and the exponent, s, of the density distribu- 
tion as well as the refractive index of the grain material 
m, we find that the measured polarization, colours and 
brightness distribution in the disc are best reproduced by 
a model in which the grains are larger than interstellar 
grains (the minimum grain size a m ; n = 0.15 /zm). The 
value of the maximum size is ill determined because it 
has little influence and was taken to be a max = 100 /im. 
The best-fit exponents of the power laws are q — 3.2 and 
s = 3. The grains have in the R-band a refractive index 
ran = 1.152 — 0.005i. Such a value is roughly appropriate 
for porous grains where half of the volume is ice and the 
other half vacuum, or where 24% of the volume consists 
of silicate and the remaining 76% of vacuum. 

Key words: Stars individual: j3 Pictoris - circumstellar 
matter - planetary systems - ISM: dust - polarization 



1. Introduction 

In the search for extrasolar planetary systems, f3 Pictoris 
and other Vega-type stars have become very popular and 
our present knowledge is summarized in several reviews 



(Backman & Paresce |1993|; Vi dal-Madjar & Ferlet [1994 ; 
Artymowicz 1994, 1996, 1997). In the last review, Arty- 



mowicz (1997) analyzed three principal components of the 



(3 Pictoris system: the star itself, the circumstellar dust 
and the circumstellar gas. Although there is an improved 
determination of the distance towa rds (3 Pic from the Hip- 
parcos satellite (Crifo et al. 1997 ), the main conclusions 
of the previous analyses remain in force. 



In general, information on dust comes from measure- 
ments of extinction, polarization, scattered light and ther- 
mal emission. For scattering, the geometrical relation be- 
tween source, scatterer and observer is essential. Whereas 
it is ill determined in reflection nebulae and allows only 
a very rough derivation of the grain properties, the geo- 
metrical configuration in the circumstellar disc of (3 Pic is 
clear and simple. In the case of (3 Pic one observes: 



i) no circumstellar extinction (Crifo et al. 1997[ ); 
ii) very weak polarization of the star itself (Tinbcr- 
gen [L979D; 



Hi) scattered light from the nearly edge-on disc (Smith & 
Terrile |1984t Kalas & Jewitt |1995| ; Mouill et et al. [1997D 
and its polarization (Gledhill et al. 1991; Wolstencroft 
199§ ; 



iv) 



et al. 

an infrared excess (Aumann et al. 1984) which extends 
up to 1300/im (Chini et al. 1991 ); there are al so mi d 
IR images (Lagage & Pantin 1994; Pantin et al. 1997). 



Many numerical calculations have been performed to 
explain the IR emission of (3 Pic (see detailed discussion 
in Li & Greenberg 1998). Depending on the wavelength 
range considered, the observations were reproduced using 
compact or fluffy grains with sizes ranging from smaller 
1 /im to up to 1 mm. 

With respect to modelling the scattering and polar- 
ization of light from the disc of f3 Pic, Artymowicz et 
al. qi98SD, Kalas & Jewitt Ql995[ |1996|) and Pantin et 



al. (1997) considered scattering only at one wavelength. 
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They assumed either isotropic or anisotropic scattering 
without computing the asymmetry factor g. Scarrott et 
al. ( 1992| ), on the other hand, applied Mie theory and 
treated multi-wavelength scattering, however, the polar- 
ization was o nly ca lculated at the then available R-band. 
Artymowicz ( 1997 ) reproduced the observation in the V- 
band employing the empirical phase and polarization func- 
tion of zodiacal and cometary dust. 

In this paper, we model scattering and polarization 
at all observed wavelengths with particle cross sections 
computed from Mie theory. As a result, we are able to 
constrain the properties of the grains and to exclude cer- 
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tain choices of dust models which were hitherto thought 
possible. 



2. Observations of polarization and colours 

Imaging polari metry of the Pic disc was performed by 
Gledhill et a l. (1991 ) in the R waveband and by Wolsten- 
croft et al.( 1995|) in the B, V, R and I wavebands. The 
polarization patterns are centro-symmetric and indicate 
that (3 Pic itself is the illuminating source. The observa- 
tional data are collected in Fig. |] and suggest the asym- 
metry in the polarizing behaviour between the SW and 
NE sides, especially in the I-band. 

The wavelength dependence of polarization P(X) is 
typical of reflection nebulae and characterized by the in- 
crease of the polarization degree with wavelength. The in- 
crease is rather weak in the NE side and well pronounced 
in the SW side. 

The polarization vectors are oriented perpendicular to 
the d isc (G ledhill et al. |199l| ; Wolstencroft |199c$ . Tinber- 
gen (1979) included f3 Pic in the list of zero polarization 
standard stars. Using a large diaphragm which included 
the disc he found that the mean degree of polarization over 
the wavelength range 0.4 - 0.7 /imisP = 0.020 ± 0.008%. 
This value is compatible with the absence of material 
in front of the star under the assumption of a maxi- 
mum polarization efficiency (P/Ay — 3%) and implies 
Ay — 0. 006 m ag. The small polarization observed by Tin- 
bergen ( 1979 ) should be due to dust scattering in the disc. 

In ordinary reflection nebulae, the colour of the scat- 
tered light is usually bluer than the illuminating star 
at sm all off sets and redder at large distances (Voshchin- 
nikov 1985 ). Unfortunately, the disc colours of (3 Pic were 
observed only out to ~12" of fset (P aresce & Burrows 1987 ; 
Lecavelier des Etangs et al. 1993 ). The data of Fig. |2j in- 
dicate that the disc has the same colour as the star or is 
slightly redder. The colours do not depend on position, 
the only exception being the innermost point (tp = 2'.'5). 
It falls into the central gap which is presently the subject 
of various speculations. 

We thus conclude that the properties of the scattered 
light are not peculiar. We therefore believe that the usual 
approach to the study of reflection nebulae can also be 
applied to (3 Pic. 



3. Modelling 

3.1. Polarization mechanism and disc model 

The variation of the position angle and the degree of po- 
larization in the disc of (3 Pic speaks in favour of light 
scattering by spheres or arbitrarily oriented non-spherical 
particles. The disc is seen nearly edge-on. As it is opti- 
cally thin, we can choose a model of single scattering in 
an optically thin medium. In this case, the radiation es- 
caping from the disc at the angular distance tp from the 
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Fig. 1. Results of the multi-wavelength measurements of 
the polarization along the disc of (3 Pic. 

star is the integral along the line of sight over some range 
of scattering angles (61, 180° — ©1), where 



©1 = arcsin 



Y 



Yu 



(1) 



and </? m ax denotes the maximum angular siz e of th e disc. 
We will adopt ip max = 45"0 (Kalas & Jewitt |1995[ ). With 
growing tp, the range of scattering angles, which is always 
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Fig. 2. Observed colours along the disc of /3 Pic. 

centered at 90°, becomes narrower. A slightly tilted disc 
orientation (^10°) would not change the picture of light 
scattering significantly. 

3.2. Step 1: polarization diagrams vs scattering angle 

As a first step in the modelling, we compute the polariza- 
tion and scattered intensity of an elementary volume. In 
Mie scattering by spherical particles, such polarization at 
a given wavelength A depends on the refractive index of 
the particle, m\ — n\ — k\i, and its radius a, 



P(m,a,\,G) 



ii(m, a, A, 9) — izim, a, A, 0) 
ii(m, a, A, 0) + ii(m, a, A, 0) 



(2) 



Here i\ and ii are the dimensionless intensities which de- 
termine the radiation scattered perpendicular and parallel 



to the scattering plane, respectively (van de Hulst 1957; 
Bohren & Huffman |1983| ). 

Figure |3| shows the polarization diagrams at A = 
0.70 fim (R-band) for which the observational database 
is largest (Fig. |l|). The particles considered are a mix- 
ture of silicate and ice. They could be porous, i.e. contain 



some volume fraction of vacuum. Porous grains have been 
discussed many times for the disc of (3 Pic (see Artymow- 
icz 1997; Li & Greenberg 1998). The refractive indices 
used are specified in Table [lj. Note that we chose the vol- 
ume fractions of the grain constituents in such a way that 
several of the refractive indices in Table |l| are identical, 
although they refer to grains of different chemical com- 
position. The cross sections of grains with a silicate core 



and an ice mantle were computed from Guttler's (1952) 
theory for two-layered spheres. 

It is important to note that all polarization diagrams of 
Fig. m are similar and independent of the particle compo- 
sition and structure. Small grains belong to the Rayleigh 
domain where the polarization has the well known bell- 
like shape with the maximum at = 90°. Polarization 
diagrams for very large spheres are also simple. They re- 
semble smooth curves with the maximum at ~ 70°. 
In both cases, the polarization does not change sign and 
reaches ~100%. For particles of intermediate size, polar- 
ization reverses sign repeatedly and a ripple-like structure 
is seen. 

This behaviour of the curves -P(A) in Fig. || reflects 
the general principles of light scattering by small par- 
ticles and does not principally change for non-spherical 
bodies, like spheroids, cylinders, bispheres, or fluffy aggre- 



gates arbitrarily aligned in space (Mishchcnko ct al. 1996 



Kolokolova et al. 1997 ); Lumme et al. 1997 ). Moreover, 
Lumme et al. ( 1997 ) demonstrated that Mie theory of 
spheres reproduces the polarization diagrams for complex 
particles rather well, except sometimes for forward and 
backward scattering. 

The polarization measurements in f3 Pic were made 
at offsets from ip = 8'.'5 to 31'.'2 (see Fig. [l]). According 
to Eq. (0), the scattering angles vary in the limits from 
(11°-169°) to (44°-136°). With these limits, we conclude 
with the aid of Fig. ^[ the polarimetric observations of (3 
Pic cannot be explained by light scattering of very small 
(< 0.1/im) or very large (> 10/im) particles only because 
they produce the polarization too high for the scattering 
angle range. 

3.3. Step 2: polarization diagrams vs particle size 

As the second step of the modelling procedure we include 
averaging along the line of sight (scattering angles) at fixed 
angular distance ip. We adopt that the number density of 
dust grains in the disc has a power-law distribution 



n d (r) 



0, r < r , 



(3) 



Here, ro is the radius of the central hole where the density 
is equal to no and r out is the outer radius of the disc. The 
latter is given by r ou t = Dip max = 19.28 • 45 » 87 AU, 
assuming a distance D — 19.28 pc (Crifo et al. 1997 ). The 
central cavity reaches out to ip <, tp w 6'.'0 (Artymowicz et 
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Table 1. Grain refractive indices ran at A = 0.70 



Refractive index Particle composition 



1.508 - 0.020i 
1.310 - O.OlOi 
1.152 - 0.005i 



1.715 — 0.030i compact silicate grains 

composite grains: mixture of 50% silicate + 50% ice 
porous grains: mixture of 72% silicate + 28% vacuum 

compact grains of dirty ice 

porous grains: mixture of 45% silicate + 55% vacuum 

porous grains: mixture of 50% ice + 50% vacuum 
porous grains: mixture of 24% silicate + 76% vacuum 
1.715 — 0.030i/1.310 — O.OlOi core-mantle grains: silicate core + ice mantle 



Optical constants for silicate are from Draine (1985), for dirty ice from Greenberg (1968). 

The optical constants for composite and porous grains were obtained from the Bruggeman rule (Bohren & Huffman 1983| ). 



al. |1989D .F1 With these numbers, r « 120 AU and the ratio 
of the inner to outer radius equals </?o/Vmax = fo/fout ~ 
0.13. 

The polarization degree can be expressed as 



P(m, a, A, ip, s) 



K 

<Ij>= ~s+l 



<h> - <h> 
<h> + <h> ' 



(4) 



7T-01 



ij (m, a, A, 0) sin s 6 d9 , (5) 



where the intensities Ij (j = 1, 2) depend on m, a, A, tp, s, 
and K. is a constant. If ip < 6'.'0 and sinG > • y . 

' 0.13 tpmax 

the line of sight intersects the hole. 

Examples of polarization diagrams are shown in Fig. |] 
at the offset ip — 10'.'5 for different types of particles and 
various values of s. The polarization degree observed in 
the R-band at this offset angle in the NE and SW exten- 
sions is shown by horizontal dashed lines. From this figure 
we again conclude that grains of intermediate size have to 
be included into consideration to explain the observed po- 
larization. Interestingly, the polarization observed at given 
offset may be explained for any density distribution of par- 
ticles: as evident from Fig. [|a (the parameter s determines 
only the polarization level at low and large size limits). 

The curves plotted in Fig. ||b look even more similar if 
we use the product \m — l\a as an argument. That this is 
so is well known for the extinction efficiency (e.g. Green- 
berg 1968). The corresponding polarization diagrams are 
shown in Fig. [|. Thus, considering light scattering in the 
continuum, the effects of refractive index and particle size 
cannot be separated. In other words, we can only estimate 
\m — l\a, the product of refractive index times particle size, 
from observations at one wavelength. 



1 The adaptive optics observations made by Mouillet et 
al. (1997) show that the scattered light is present at angu- 



lar distances up to ip « 1. 5. Perhaps, this is the light from 
the outer disc scattered in almost forward and backward di- 
rections. However, because the polarization observations are 



In more realistic models, we must also average over the 
size distribution n<i(a). Afterwards the spikes in polariza- 
tion seen in Fig. || disappear. We use the power-law 



rid (a) 



(6) 



with minimum and maximum radii a m ; n and a max , re- 
spectively. Figures || and Q are analogous to Fig. [| and 
H but include an average over the size distribution. The 
upper radius and the exponent are kept fixed (a max = 
100 /im, q = 3.5) and a m i n is being varied. 

The changes of a m - ln are most important with respect 
to polarization. The increase or decrease of the maximum 
size has no influence if m — l|a max 10. So if very large 
particles are present in the disc of f3 Pic, they are not 
visible as scatterers at waveband R. Changes in q do not 
affect the overall picture and will be discussed in Section 
3.6. Figures || and ^ contain plots for all refractive indices 
listed in Table 0. The fi gures show that in order to ex- 
plain the observed polarization one needs to keep a m i n in 
the range 0.6 /im Si a m i n 10 /im. However, very porous 
grains with small values of a m i n may, in principle, be used 
to model the observations as well. 

3.4- Step 3: polarization vs angular distance from star 

In the next step, we model the angular distribution of 
polarization using the minimum cut-off in the size dis- 
tribution as estimated from Fig. ^. The curves in Fig. || 
demonstrate that the same observations may be satisfacto- 
rily fit with particles of different composition by changing 
only the parameters a m i n and s.[] Note that the variations 
of s and a m ; n influence the polarization mainly at small 

made outside of this radius, this does not impact on the result, 
given in the single scattering case. 

2 Actually, the choice of the parameter s is determined by 
the radial brightness distribution of the scattered light, I(tp). 
For the outer parts of the disc of (3 Pic Up > 6'.'0), I(<fi) oc tp 1 
with 7 = 3.6 - 4.3 (Kalas & Jewitt |l995|). The values of 7 for 



the models considered are given in Table 
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Fig. 3. The degree of linear polarization for various kinds 
of particles. Their composition and refractive indices are 
given in Table [|. The peak of the curve for a — 0.1 is 
always slightly shifted to larger values of scattering angle 
relative to those for a = 0.01 /im. 
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Fig. 4. (a) Upper panel: polarization vs particle size at 
the angular distance (p = 10. 5 for various exponents, s, of 
the dust density distribution. The grains are compact sil- 
icates with ttir as indicated. The horizontal dashed lines 
show the observed polarization Pr in the NE and SW 
sides, (b) Lower panel: as (a), but now for a fixed den- 
sity distribution (s = 2) and particles of various refractive 
indices. 
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Fig. 5. As Fig. 4b, but with polarization plotted vs the 
parameter \mn — l|a. 

and large offset angles. (We adopted a max = 100 /xm and 
q = 3.5.) Therefore, the polarimetric observations of the 
extended sources (nebulae, circumstellar shells, galaxies) 
at one wavelength do not allow to determine the chemical 
composition of particles. Numerous other sets of refractive 
indices and size parameters may also explain the observa- 
tions. As an example, let us consider the model of Scarrott 
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Fig. 6. Polarization vs a m i n at the angular distance <p = 
10'.' 5. The upper size limit equals a max = 100 /im and the 
exponent of the size distribution q = 3.5. 
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et al. (1992). They reproduced the angular dependence of 
the scattered light and polarization in the R-band ob- 



served by Glcdhill ct al. (1991) with silicate particles of 
mfl = 1.65 — 0.05z and a m ; n = 0.01 /im, a max = 3/um, 
q = 4.0, s = 2.75. However, the disc colours in their model 
are too red compared to those observed (from — 0™72 in 
(B-V)*_ dis c to -l m 24 in (V-I)*-«Ubo). 

3.5. Step 4-' polarization vs wavelength 

With the parameters found from modelling of the polariza- 
tion in the R-band (see Fig. ||), we calculated the wave- 
length dependence of polarization at the offset distance 
(p = 10'.' 5. We also checked the colour excesses given by 
the corresponding model (see Table |^) and compared them 
with observations. 

The results are shown in Fig. ^ and Table |2|. It is clear 
that despite our thorough procedure and the successful 
fits of the changes in polarization and brightness in the 
R-band with angle <p (see Fig. || and Table ||), the de- 
pendence P(A) cannot be well explained for any refractive 
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Fig. 8. Polarization in the disc of (3 Pic as a function 
of the angular offset for the R-band. Solid lines present 
models with different refractive indices (see Table [[]). In 
all calculations, a max = 100 fim and q — 3.5. The other 
model parameters (a m i n and s) are shown in the panels. 



index. From Table ^ it also follows that the colour index 
V-I is crucial for the choice of the model. The colours in 
the models 1-4 are too red and these models must be 
rejected. In principle, the correct polarization behaviour 
(i.e. its growth with wavelength) can be reproduced for all 
refractive indices of Table |l| if the minimum size a m ; n is 
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Table 2. Observational and theoretical colour excesses in the disc of (3 Pic at <p = 10'.'5. 
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Colour 


excess 


Observations 


Model 1 


Model 2 


Model 3 


Model 4 Model 5 


(B-V)* 


— disc 


-CP21 ± 0P2Q 


-0 m 217 


-CT213 


-0^221 


-0^228 m 034 


(V-R)* 


— disc 


-(ri7±0 m 20 


-(T200 


-(F226 


-(F228 


-(T252 OP061 


(V-I)*- 


disc 


m 01 ± m 20 


-(T313 


-(T344 


-0 m 354 


-CT400 (T105 


Model 1: 


m R = 


1.715 - 0.030i, a m 


n = 0.7 fxm, 


S = 2 ( 7 


= 3.4, A = 


0.67, g = 0.80). 


Model 2: 


mn = 


1.508 - 0.020i, a m 


n = 1.2 jum, 


s = 2 ( 7 


= 3.4, A = 


0.66, g = 0.86). 


Model 3: 


m R = 


1.310 - O.OlOi, a m 


n = 2.4 ^m, 


s = 3 (7 


= 4.2, A = 


0.67, g = 0.91). 


Model 4: 


m fl = 


1.152 - 0.005i, a m 


n = 5.5 fJXCl, 


s = 3 (7 


= 4.4, A = 


0.67, g = 0.96). 


Model 5: 


mfl = 


1.152 - 0.005i, a m 


n = 0.2 /urn, 


s = 3 (7 


= 4.4, A = 


0.89, g = 0.82). 



All models have a max = 100 fim and g = 3.5. 

7 - exponent in the brightness distribution, A - single scattering albedo, g - asymmetry factor. 



lowered. But even if we were to fit the dependence P(X), 
very small particles would yield unacceptably blue colours. 

So in order to explain all colour and polarization ob- 
servations we need to reduce the lower limit of the size 
distribution, a m j n , and the refractive index m. 

3.6. Step 5: polarization and colours: whole picture 

Finally, we suggest a fit to the angular dependence of 
the polarization (in all wavebands) and the colour ex- 
cesses (see Figs. |l0| and [ll]). In order to improve the fits, 
we varied the parameters a m i n and q around the values 
given in Table || for model 5. Unfortunately, there re- 
mains some problem for the position with the minimum 
distance ip = 2'.'5. The colour index B-V can be mod- 
eled only if the grains have very narrow size distribution 
around a ~ 15 /im, but with such grains the polarization 
at observed offsets becomes too large (see Fig. ^). 

4. Concluding remarks 

All models presented here were constructed with a small 
number of parameters which are rather well determined. 
Our favourite model presented in Figs. [HJ ^ has param- 
eters: 

m R = 1.152 - 0.005i, 

Amin = 0.15 (MID-, 
Omax = 100 /Ltm, 

9 = 3.2, 
s = 3. 

The model fits reasonably well the observed angular de- 
pendence of polarization (in four wavebands), three colour 
excesses and the radial brightness distribution. The calcu- 
lated values of the exponent in the brightness distribution, 
single scattering albedo and asymmetry factor in the It- 
band are: 

7 = 4.4, 

A = 0.88, 

g = 0.79. 

Still better fits could be obtained assuming different dust 



properties on either side of the disc. Previous estimates of 
the grain albedo in the disc of [3 Pic based on visiblc/IR 
models gave A = 0.5 ± 0.2 (Artymowicz 1997), which is 
smaller than our value for very porous grains. But Arty- 
mowicz et al. ( [198SD also derived A ~ 0.9. 

As usual, the albedo A of the grains was estimated in- 
dependently of the scattering phase function. Evidently, 
the latter cannot be isotropic in the visual because so 
far all suggested grain sizes imply forward scattering. The 
value of our model (g w 0.8) is larger than previous ones. 
Using t he H eyney-Greenstein phase function, Kalas & 
Jewitt ( 1995 ) found that the brightness distribution ob- 
served in the disc of [3 Pic in the R-band may be repro- 
duced with 0.3 < \g\ < 0.5. Note t hat m i rror p articles 
(with g < 0) used by Kalas & Jewitt (|l995| , |l99rj ) are as- 
tronomical nonsense because no grain material (with the 
exception of very small iron particles) gives at visual wave- 
lengths a negative asymmetry parameter. Even pure con- 
ductors (m = oo) have g w —0.2 only (van de Hulst 1957 ). 



10L [Ll]) we only spec- 
It was taken arbi- 



In our favourite model (see Figs, 
ify the refractive index (see Table |) 
trarily and is not motivated from the physical point of 
view. From other side, the porous particles may be easily 
kept in the shell because the radiation pressure acting on 
them is smaller: we found that for silicate grains of radii 
a <; 2.0 /im with 76% porosity (Table Q) the radiation 
pressure force is smaller than gravitational force. 

Very porous aggregates were adopted by Li & Green- 



berg (1998) to explain the IR observations. Mie theory 



can, of course, not be applied in such cases. Although Li & 
Greenberg used a standard mixing rule to treat the fluffy 
aggregates, the underlying theory has not been developed 
for particles with sizes larger than the wavelength nor for 
computing the scattering of light. 

We modeled scattering and polarization in the disc of 
(3 Pic at all wavelengths where it has been observed. Be- 
cause the properties of the scattered light are not pecu- 
liar, we could compute particle cross sections from Mie 
theory. Our models exclude the possibility that the grains 
are compact spheres, or that they are all very small or 
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Fig. 9. The wavelength dependence of polarization in the 
disc of (3 Pic at the angular distance ip = lf/.'5. The ob- 
servations in the SW and NE sides and their errors are 
plotted by dashed lines. The parameters of the model are 
the same as in Fig. 8. 
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Fig. 10. The polarization observed in the disc of (3 Pic. 
The solid lines correspond to the theoretical models for 
which inn = 1-152 — 0.005i, a m i n = 0.15 ^m, a max = 
100 /itm, q = 3.2 and s = 3. 



all very large grains. Instead, the particles must be rather 
porous (degree of porosity ^ 50 %) and have a lower cut- 
off in the size distribution as small as a m { n ^0.2 ^m. This 
limit is larger than the mean size of interstellar grains. Our 
model can be further verified by computing the resulting 
IR fluxes and check how well they agree with observations. 
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Fig. 11. The colour excesses observed in the disc of (3 Pic. 
The theoretical model has the same parameters as in 
Fig. 10. 
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